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We consider S — 1/2 antiferromagnetic Heisenberg chains with alternating bonds and 



quenched disorder, which represents a theoretical model of the compound CuCl2a;Br; 



Xl-x) 



(7- 



pic)2. Using a numerical implementation of the strong disorder renormalization group method 
we study the low-energy properties of the system as a function of the concentration, x, and 
the type of correlations in the disorder. For perfect correlation of disorder the system is in the 
random dimer (Griffiths) phase having a concentration dependent dynamical exponent. For 
weak or vanishing disorder correlations the system is in the random singlet phase, in which 
the dynamical exponent is formally infinity. We discuss consequences of our results for the 
experimentally measured low-temperature susceptibility of CuCl2a;Br2(i_i) (7 — pic)2. 
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1. Introduction 

Low-dimensional quantum systems (antiferromagnetic 
spin chains, ladders, two-dimensional systems, etc.) are 
fascinating objects which have been investigated inten- 
sively in experimental and theoretical works. From the 
theoretical point of view these systems exhibit several 
unusual properties, for example quasi-long-range order, 
topological string order, a spin liquid state, quantum 
phase transitions, etc.. Many of these features can be 
illustrated by the 5 = 1/2 antiferromagnetic Heisenberg 
chain, which is defined by the Hamiltonian: 



i+1 1 



(1) 



in terms of the spin-1/2 operators, Si, at site i. The ho- 
mogeneous model with Ji = J is gapless and there is 
quasi-long-range order in the ground state, i.e. spin-spin 
correlations decay algebraically,^ {SlSl_^_^) ~ r~^. Here 
T — x,y,z and (...) stands for the ground-state expecta- 
tion value. Introducing bond alternation such that Ji — J 
at i = even and Ji = J, (ad > 0) at i = odd, the sys- 
tem with Qfrf 7^ 1 is in the dimerized phase, ^ in which 
spin-spin correlations are short-ranged and the excita- 
tion spectrum has a finite gap. In the bond alternating 
model the quantity 5d — Ina^ serves as a quantum con- 
trol parameter and the quantum critical point is located 
at 51 = 0. 

Quenched (i.e. time independent) disorder has a pro- 
found effect on the low-temperature/low-energy proper- 
ties of quantum systems,"^ both at the quantum criti- 
cal point and in an extended region of the off-critical 
regime-the so called quantum Griffiths phase.*' ^ In the 
disordered version of the uniform S — 1/2 antiferro- 
magnetic Heisenberg chain defined in eq.(l), the cou- 
plings, Ji, are independent and identically distributed 
random variables. A detailed study by Fisher^ using 



an asymptotically exact strong disorder renormalization 
group (SDRG) method^ revealed that the ground state of 
the random model is the so called random singlet phase, 
in which singlets are formed between spin pairs which 
could be arbitrarily far apart. Average spin-spin corre- 
lations, which are dominated by rare regions, are quasi- 
long-ranged: [(S'iS'i+r)]av ^ 1 where [. . . ]av stands for 
averaging over quenched disorder. On the other hand, 
typical correlations are much weaker and decay asymp- 
totically as ln((S'iS'i+r)typ) ^ r^/"^. In addition, the dy- 
namical scaling in the random singlet phase is anoma- 
lous. The length scale L and the energy scale A, mea- 
sured by the lowest gap, are related logarithmically: 



InA-L"^, -0 = 1/2. 



(2) 



In the random bond alternating Heisenberg model bonds 
at even (Jc) and odd (Jo) sites are taken from different 
distributions and the quantum control parameter is de- 
fined as 



5 = [In Jc]av - [In Jo] a 



(3) 



For 6^0 the model is in the random dimer phase, in 
which spatial spin-spin correlations are short ranged.* 
Within a range of finite dimerization \S\ < So > 0, dy- 
namical correlations are however still quasi-long-ranged, 
which is related to the fact that the system is gapless. 
This region is called Griffiths phase. In a finite chain of 
length L the typical gap scales asymptotically as: 

A - . (4) 

Here the dynamical exponent, z, is a continuously vary- 
ing function of the control parameter, S. On the border 
of the Griffiths phase, 6 — da, the dynamical exponent 
vanishes, whereas close to the random singlet phase it 
diverges with 6 as:^ 



(5) 
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In the random singlet phase as well as in the Grif- 
fiths phase, thermodynamic quantities such as the low- 
temperature susceptibility, xi^), the low-temperature 
specific heat, C{T), and the low-field magnetization, 
M{h), at zero temperature are singular. For example the 
susceptibility in the random singlet phase behaves as: 

^^^^ ~ T^r? ' ^""^ 

whereas in the Griffiths phase it is given by: 

X(T)~T-1+^ p = l/z. (7) 

Recent experimental stiidics on the compound 
CuCl2a;Br2(i_j,)(7 — pic)2 show another type of competi- 
tion between bond alternation and randomness.^' Here 
the .X = compound, CuBr2(7 — pic)2. is a homoge- 
neous 8 = 1/2 antiferromagnetic Heisenberg chain in 
which the Gu — Gu bond is bi-bridged by two Br atoms: 
Cu < > Gu and the coupling constant is given by 
J" = 20. 3K.^ The x = 1 compound, GuGl2(7 — pic)2, is a 
bond alternating 5 = 1/2 antiferromagnetic Heisenberg 
chain in which two kinds of Gu < q[ > Gu bonds alter- 
nate along the chain. The bond alternation is induced 
by a freezing transition of the rotational motion of the 
methyl-group at bQK .^"^ The experimentally measured 
coupling strengths are J = 13. 2K and Ja, with a = 0.6. 
In the mixed compound with a small finite concentra- 
tion of Br [x 7^ 1), atoms connected with alternating 
Gu < ^[ > Gu bonds form a cluster and different clus- 
ters are separated by bonds with Br — Br and/or Br — Gl 
bridges. The strength of the Gu < > Gu bond has 
been estimated from the theoretically calculated mag- 
netization curve as J' = 1.3J.^^ In ref. 9 the data for 
the low-tempcrat\ire susceptibility of the diluted system 
show an algebraic temperature dependence, x ^ T'^^^, 
which is compatible with the expected behavior in the 
Griffiths phase as given by cq.(7). For a wide range of 
the concentration, a;, the measured effective exponent is 
about /?exp ~ 0.5 — 0.67, and shows only a weak concen- 
tration dependence. 

2. Theoretical Model 

A theoretical model for GuGl22;Br2(i_j;)(7 — pic)2 was 
presented and numerically studied in ref. 12. It is not 
obvious from the information one can extract from the 
experiments whether the rotational order of the methyl- 
group remains long-ranged in the presence of the dilu- 
tion by Gu < ^1 > Cu or Gu < > Gu. There- 
fore two models were introduced in:"'^^ i) The fixed parity 
model in which the rotational order is assumed to be 
perfectly long ranged, i.e. the J (a J) bonds stay in the 
same parity position, say at i = even (i = odd), in any 
bond-alternating cluster. Introducing a parity parame- 
ter. Pi, given by Pi = 1(— 1) if the Gu < ^{ > Gu bond 
has a value of J (aJ), one can describe the parity cor- 
relations in the fixed parity model as PiPi+21 = 1 and 
PiPi+2i-i = —1; fi) The random parity model in which 
the rotational long-range order between two clusters of 
alternating bonds is assumed to be completely destroyed 
by the dilution. In this case, if i and i + 21 {i + 21 — 1) 
refer to different clusters, the average parity correlations 
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vanish: \piPi+2i]av = \PiPi+2i-i]a.v = 0. 

In this paper we consider a more general model 
in which rotational long-range order is partially de- 
stroyed by dilution, so that correlations between pari- 
ties in two different clusters are given asymptotically as 
[PiPi+2i]av = -[PzPi+2i-i]e.v ~ Here we recover the 
fixed parity and the random parity models in the lim- 
its p = and p ^ 00, respectively. We note that some 
aspects of the effect of correlated disorder on quantum 
systems is studied in ref. 13. Here we summarize the val- 
ues of the coupling constant in eq.(l) in the following 
way: 

( Ja^^-P'y^ with prob. x'^ 
Ji= I J' with prob. 2x{x - 1) . (8) 

[ J" with prob. {x - if 

The parity correlations within a cluster of alternating 
bonds are given by PiPi+j = (—!)■' and the parity corre- 
lations for bonds in different clusters are defined above 
for the different models. 

The theoretical results for the low-energy behavior of 
the fixed parity and the random parity models obtained 
from the density matrix renormalization group (DMRG) 
method^^ are not fully consistent with the measured low- 
temperature properties of the GuGl2rEBr2(i_a;)(7 — pic)2.^ 
In ref. 12 the random parity model is found to be in 
the random singlet phase, in which the susceptibility 
exponent is /? = and is much lower than the exper- 
imental value /3exp ~ 0.5 — 0.67. For the fixed parity 
model the susceptibility exponent is found to be too 
large: /? > 1. Another problem with the latter model 
is that the DMRG analysis could not be performed for 
X < 0.4 due to strong finite size effects. 

In the present paper we revisit the models of the com- 
pound GuGl2a:Br2(i_a;)(7 — pic)2. Our study is different 
from ref. 12 in two respects. Firstly, we consider a more 
general model in which the effect of disorder correlations 
are taken into account. Secondly, we use a numerical im- 
plementation of the SDRG method. This method usually 
gives very accurate results for the form of singularities, 
in particular for the dynamical exponent. This method 
has been successfully applied to clarify the low-energy 
singularities of more complicated systems, such as ran- 
dom spin ladders^^ and two- and three-dimensional ran- 
dom Heisenberg antiferromagnets.^^ The major advan- 
tage of this method is that one can consider large systems 
L ~ 1000 — 2000 with good statistics compared with the 
DMRG method. 

3. Numerical results 

3.1 Numerical method 

The SDRG method for random AF spin chains pro- 
ceeds as follows: During the renormalization process, one 
first identify the strongest bond, say J23 = O, which con- 
nec;ts sites 2 and 3. If the neighboring bonds are much 
weaker, J23 ^ J12, J34, then the spins 2 and 3 form an 
effective singlet and can be decimated out. The new cou- 
pling between the remaining sites, 1 and 4, is obtained 
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within a second-order perturbation calculation as: 

j = \'h^. (9) 

By repeating this decimation process, we gradually re- 
duce the energy scale, 17. At the fixed point, the energy 
scale is give hy VI. ^ \i only a small fraction 1/L of 
sites is active (i.e. not yet decimated). In our computa- 
tions we consider finite periodic chains with L = even 
sites and decimate until the last pair of sites having a < 
gap A J7. From the distribution of the logarithms of 
the gap: 



agrees very well with the known exact result^^ for \/z 
which is given by the formula [{Je/ JoY^' ''^] = 1, i.e. by 
the solution of the equation (xJ^^/^^ + (1 — a;)(aJ)^^/^)) • 
((1 - a;) + a;(aJ)(-i/^)) = 1. 



(a) 



(b) 



P(lnA) 



A^O, 



(10) 



we obtain the the dynamical exponent, z. In the random 
singlet phase, the dynamical exponent is formally infi- 
nite, as described in eq.(2), and the appropriate scaling 
combination in a finite system is given by: 



ln[L'^P(lnA)] ~ /(L^'^lnA) 



(11) 



The SDRG method as outlined above is expected to pro- 
vide asymptotically exact results both in the RS phase^ 
and in the Griffiths region. 
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Fig. 2. Distribution of the energy gap for the fixed parity model 
for two concentrations: (a) x = 0.6, (b) x = 0.8. The asymptotic 
slope of the curves, indicated by broken straight lines, gives the 
susceptibility exponent, /3. 
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Fig. 1. Test of the SDRG method on random dimerized XX- 
chains, with a = 0.6, x = 0.6. The slope of the distribution of 
the gaps in a log-log scale agrees well with the known exact result 
for the dynamical exponent, which is given by the slope of the 
broken straight line. 

In the present model the randomness is discrete and 
at the starting point of the renormalization there are 
several couplings with the same value of f2. In this de- 
generate situation we randomly choose the actual cou- 
pling to be decimated. After a sufficiently large number 
of renormalization steps we will have a (quasi)continuous 
distribution. To illustrate the correctness of the SDRG 
procedure for discrete randomness, we consider the ran- 
dom dimerized XX-c\\.am^^ in which the couplings take 
the value J or aJ and a fraction, x, of the odd (even) 
couplings are J (aJ). To be close to our model in eq.(8) 
we took a — 0.6 and x = 0.6 for the XX-chahn. In Fig. 1 
the distribution of the gaps for different finite systems is 
shown in a log-log scale, the asymptotic slope of which 



Fig. 3. Average gap of the fixed parity model with different con- 
centrations, X, as a function of the size, L. The asymptotic slope 
of the curves in log- log scale defines the dynamical exponent, 
z = l/p. Note the strong finite size corrections: for sizes used in 
DMRG (In L <b) we are not in the asymptotic regime. 



3. 2 Fixed parity model 

Next, we turn to our model defined in eq. (8) and start 
with the fixed parity case. The distribution of the gaps 
for two different values of the concentration, x — 0.6 
and 0.8, are shown in Fig. 2. Evidently, the exponent 
(3, given by the asymptotic slope of the distributions, is 
finite and concentration dependent. It can also be ex- 
tracted from the scaling behavior of the average gap via 
[InAJav — const -|- /9~^lnL. As seen in Fig. 3 one can 
obtain a reliable estimate of /3 for a; > 0.4. For smaller 
concentrations, even the largest system size L = 2048 
is not yet in the asymptotic regime. The estimates for 
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/3 are depicted in Fig. 4. They are systematically lower 
than the DMRG results in ref. 12, which is probably due ^ 
to finite-size eflFects. Using the DMRG method, ref. 12 
reports results up to InL < 5, which is, even for com- 10"^ 
paratively large concentrations, not in the asymptotic T 10"^ 
regime according to Fig. 3. 
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Fig. 4. The susceptibility exponent, /3 = 1/z, for the fixed par- 
ity model as a function of the concentration. For small x the 
expected quadratic dependence is denoted by the dotted curve. 

As can be seen in Fig. 4, for small concentrations there 
is a quadratic dependence: (3 ^ x? . which can be un- 
derstood as follows: In the fixed parity model the aver- 
age value of the log-couplings in the odd and even posi- 
tions arc different, so that the dimcr control parameter 
in eq. (3) is 5 ~ Ina ^ 0. For small x, thus for small 5 
we can use eq. (5), which is compatible with the observed 
quadratic a;-dependence of /?. 

3.3 Random parity model 

For the random parity model the gap distribution for 
an intermediate concentration, x = 0.8, is shown in 
Fig. 5(a). One observes that the distribution gets sys- 
tematically broader with increasing system size, which 
is a characteristic of the RS phase. Indeed, the logarith- 
mic scaling in eq. (11) is well satisfied with ^ = 1/2, 
as illustrated in Fig. 5(b). The same logarithmic scaling 
with if) = 1/2 holds for other values of the concentration, 
too. These results can be understood by noting that odd 
and even couplings have on average of the same strength 
in the random parity model, consequently the (dimer) 
control parameter is (5 = 0. The short-range correlations 
in the dimerized sequences, especially for large concen- 
tration, does not modify this large-scale asymptotic be- 
havior. Our results for the random parity model are in 
agreement with the DMRG calculations in ref. 12. 

3.4 Correlated random parity model 

Now we discuss the results for the correlated random 
parity model. We recall that according to scaling argu- 
ments^^' the disorder correlations modify the critical 
behavior in the RS phase, provided the decay exponent 



Fig. 5. Probability distribution of the gaps in the random parity 
model with rr = 0.8 (a). Scaling collapse for different sizes using 
the scaling form in eq.(ll) (b). 




Fig. 6. Scaling plot of the gap distribution for the correlated ran- 
dom parity model for different values of p < 1. The i/j exponent 
is p-dependent and well described by the relation in eq.(12). 
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Fig. 7. Variation of the estimated exponents V'(p) for the corre- 
lated random parity model with the correlation parameter p. The 
solid line represents the linear fit in Eq.(12). 
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is sufficiently small: p < 2/i/, where v = 2 is the corre- 
lation length exponent for uncorrelated disorder. To il- 
lustrate this relation we estimate the control parameter, 
5, in a finite chain of length, L. Introducing the notation 

a = ln(J2i/J2i-i), we have 5^ = 1/4^2 [(^fj/ etf ~ 

l/L G{r)dr L~p, where G{r) is the disorder corre- 
lator. Thus 6 ~ L~P^'^, for /9 < 1, whereas for p > 1 we 
have the non-correlated disorder result, S ~ L"^/^. We 
can conclude from this consideration that for p < 1 the 
correlated random parity model has a vanishing gap and 
is in the RS phase, however with p-dependent proper- 
ties. The scaled gap distribution for four different corre- 
lation parameters, p < 1, is shown in Fig. 6, in which the 
logarithmic scaling collapse in eq. (11) can be obtained 
with different exponents, ijj{p). A plot of the ip{p) data 
is given in Fig. 7 in which the p dependence can be well 
fitted with the approximate formula: 

3-p 



P<1 



(12) 



4 

This result is exact at p = 1, in which case the standard, 
non-correlated RS phase result, tp = 1/2 is recovered. 
For a small p the exponent approaches ip « 3/4, but at 
p = we jump to the fixed parity model. The linear de- 
pendence of ■ip{p) on p in eq.(12) is in the same form as in 
the exact expression for the XX-model^^ with correlated 
randomness. 

The low-temperature singularity of the susceptibility 
of the correlated random parity model is given by the 
formula: 

in which tjj is p-dependent, as given in Eq.(12). 

4. Discussion 

We close our paper by comparing the experimentally 
measured low-temperature susceptibility of the com- 
pound CuCl2xBr2(i_a;)(7 — pic)2 with the results of our 
theoretical calculations. Wc recall that the measured sus- 



ceptibility exponent is finite, /3e: 



0.5 - 0.67, and has 



only a weak concentration dependence. These properties 
are partially compatible with the results for the fixed 
parity model, which is found to be in the Griffiths phase, 
where f3 is finite for all values of the concentration. Ac- 
cording to Fig. 4 there is a range of concentration around 
X fs 0.8 in which our results for /3 agree with the ex- 
perimental values. However, this range is rather narrow 
and our values for /3, which have a substantial concentra- 
tion dependence, are generally significantly smaller than 
/3exp- For models with non-strictly correlated parities, in 
which case /3 is found formally zero, the agreement with 
the experiment is even less satisfactory. One possible ex- 
planation of this discrepancy between experiment and 
theory is that the experimentally measured /3oxp sltc ef- 
fective, temperature dependent values, which should ap- 
proach the true behavior as T — > 0. Indeed, as seen in 
Fig. 3 the local slopes of the curves that show In A vs. 
InL, which define the effective exponent z{L) = \/(3{L), 
have a strong size dependence. This can be converted 
into a temperature dependence through O ~ T ~ . 



The corrections to /3(i) are particularly strong for small 
values of x. For moderately large sizes, InL « 5, the 
effective exponents have only a weak concentration de- 
pendence. The leading finite size corrections are of the 
form, 13[L) ~ (} + a/L, which are compatible with a tem- 
perature correction as: 



/3(r) ~ /3 + cT^ + 



(3>0. 



(14) 



In the random singlet phase with (3 = 0, both for corre- 
lated and non-correlated parity, the effective exponents 
have a logarithmic temperature dependence: 

Cl 



llnTI 



+ .... p = o. 



(15) 



which can be obtained by analyzing eqs.(6) and (13). 

Indeed, these finite temperature corrections are strong 
for the fixed parity model, particularly for small (3, i.e. 
for a small concentration, x. 



K 10 




Fig. 8. Temperature dependent length-scale, 1{T), which is ob- 
tained by performing the renormalization down to an energy 
Q = T, for the fixed parity and the random parity (inset) mod- 
els. In the log-log plot the local slope of the curves gives the 
temperature dependent effective exponent, /3(T). The slope of 
the straight line indicates the effective exponent at T = 2K. 



We tried to estimate the effective exponent, (3{T), at 
the temperature of the experimental measurement T = 
2K. For this we performed the renormalization transfor- 
mation down to an energy-scale, Q = T, and measured 
the fraction of non-decimated sites: n{fl) — l//(0). Here 
the length-scale is given by l{fl) ~ fi"^. Our results for I 
as a function of the temperature for the fixed parity and 
the random parity model are presented in Fig. 8 for dif- 
ferent values of the concentration. In the log-log plot the 
local slope of the curves is just the effective exponent, 
P{T), a.t T = Q. As can be seen in this figure the effec- 
tive exponent, /9(T), approaches its limiting value only 
if the temperature is sufficiently low. At the tempera- 
ture of the measurement, T — 2K, the asymptotic region 
seems to be still quite far. For the fixed parity model 
the effective exponent is about P = 0.4, which is prac- 
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tically independent of the concentration. This result is 
consistent with the experimental results of ref. 9. For the 
random parity model, as shown in the inset, the effec- 
tive exponent continuously vary with the temperature. 
Its value at T = 2K is about /? = 0.36, which is also con- 
sistent with the experiments. Therefore, using the avail- 
able experimental data it is not possible to distinguish 
the type of parity correlations present in the compound, 
CuCl2ccBr2(i_a;)(7 — pic)2 and the question, which type 
of model should be used for its theoretical description 
remains still open. Experimental measurements at lower 
temperatures are needed to clarify this point. 
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